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Abstract 

In this work a local projection stabilization method is proposed to solve a fictitious domain problem. The 
method adds a suitable fluctuation term to the formulation thus rendering the natural space for the Lagrange 
multiplier stable. Stability and convergence are proved and these results are illustrated by a numerical 
experiment. 
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1. Introduction 

The numerical solution of problems on smooth domains, but with complicated geometries, can be faced 
using different approaches, e.g., isoparametric elements, approximating the curved boundary by a polygonal 
one, etc. The problem becomes particularly complicated in the case the domain moves or changes shape, thus 
forcing a constant remeshing if the shape is to be tackled in time. To avoid this, and other complications, a 
fictitious domain method was proposed in [7] and analyzed in [6]. The ficticious domain approach replaces 
the original problem by a mixed one on a larger (and simpler) domain that seeks for the original variable and 
a Lagrange multiplier on the physical boundary. In the analysis given in [5] it is proved that the combination 
of piecewise linear functions for the primal variable and piecewise constants for the multiplier are inf-sup 
stable and convergent under the geometrical restriction that the mesh on the physical boundary is coarser 
than the mesh induced by the triangulation of the larger domain. This is a limitation especially considering 
that the aforementioned intersection is needed to assemble the matrix associated to the discrete problem. 
Since then, some attempts have been made to overcome this restriction, such as cut elements (cf. [SJH]), or 
XFEM approaches (cf. [9J|8]). 

In this work we propose a simple solution to this problem by means of a LPS-like stabilized finite element 
method. The starting point is the observation that the mesh for the larger domain induces a partition of the 
physical boundary. The Lagrange multiplier space built from this partition contains a subspace such that 
the combination is inf-sup stable. Then, the present approach adds a fluctuation term to the formulation 
penalizing the distance between this natural finite element space and the underlying stable pair. The analysis 
of the new method follows then an approach related to the ones treating minimal stabilization frameworks, 
such as pQ and [2J. 

The rest of the paper is organized as follows. Section [2j introduces the notations and the problem of 
interest. Then, the method is presented and its stability is proved in Section[3] Section [4] contains the error 
analysis which is illustrated by means of a numerical experiment in Section (5J 

2. Notations 

We consider ui C K 2 an open bounded domain with a Lipschitz continuous boundary 7 and outward 
normal vector n. To avoid technical difficulties we will suppose that 7 is polygonal and then it's the union 
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of ./V straight segments 71,72, ...,7jv- The analysis can nevertheless be extended with minor modifications 
to the general case. For D C M 2 , the inner product on L 2 (D) (or L 2 (D) 2 ) will be denoted by (■, -)r>- We 
adopt the usual notations for Sobolev spaces. In particular, H 2 (7) will be the space of traces of functions of 
on 7, with dual ff _ 5 (7). The duality product on £f _ 5 (7) x ff 2 (7) will be denoted by (•, -) 7 . Also, 
for 8 G [0, |] the following space will be useful in the sequel 

nf =1 ^(7i) := U € i 2 (7) : £l 7j € Jf 5 ( 7i )} . 
The problem of interest reads as follows: 

— Au = f in u> , u — g on 7 , (1) 

where / € L 2 (uj) and g € if 5 (7). The fictitious domain approach relies on the introduction of a larger (and 
simpler) domain n D w, an extension / of / to f2, and the solution of the following mixed problem: Find 
(u,X) G W := ffo(^) x #~%) s^c/i t/iat 

(Vtt,V?;)n-<A ! u) 7 + <M,u) 7 = (/,?;)n + (M,3} 7 V(v,/j)gW. (2) 

Problems ([!]) and ([2]) are linked by the fact that if (u, A) satisfies ([2|, then u| w satisfies ([I]) and A coincides 
with the jump of the normal derivative of u on 7 (see [6] for details). 

To solve this weak problem, we introduce 7^, a regular triangulation of f2 built using triangles K with 
diameter /i^, and h := maxxeTi, ft-if ■ Let 7^ be the partition of 7 induced by 7/j, this is, the collection of 
edges e such that their end points are the intersections of 7 with the edges of the triangulation Th, plus 
the angular points of 7. Let also j h be a partition of 7, whose vertices are also vertices of 7/1, with edges e 
satisfying the following (cf. [6]): there exists C > such that 3/i < |e| < Ch, for all e G 7^. Using the mesh 
regularity of 7h it is easy to see that for all e G 7^, card{e G jh ■ e C e} < C, where C > is independent 
of e and /i. 

Over these partitions we define the following finite clement spaces: 

V h := {v h G C°(fi) n J#(n) : G Pi(ff) , VK G 77,} , 

A h := {q h G L 2 (7) : qt\e G P (e) , Ve G 7/l } , 

A s := {q- h G L 2 ( 7 ) : q h \s € Po® , Ve G Tfl } , 

and := V/j x A^. The pair Vh x A^ is not inf-sup stable, while, thanks to the hypothesis on Th and 7 r, 
the pair V/, x A^ satisfies a discrete inf-sup condition (cf. j6]). 



3. The stabilized formulation and its stability 

To avoid the need to use the space A^, in this work we propose the following alternative discrete problem: 
Find (uh,Xh) G Wh such that: 

B[(uh,Xh),(vh,^h)] = (f,Vh)n + (g,fJ l h)-y V (vh, nil) G Wh , (3) 

where 

B[{u h , Xh), {v h ,pL h )] = {Vu h ,Vv h )ii - (X h ,v h } 1 + (fi h ,u h )^ + ^2 C s\e\(Xh - PX h ,^h ~ P^h)e , (4) 

C s > 0, and P : L 2 {e) — > Po(e) stands for the orthogonal projection in L 2 (e), i.e., P£|g := |e| _1 (^, l)g. 
Before heading to stability, we state the following preliminary result. 

Lemma 1. There exists f) > such that, for all fj,h G A/,. 

1 

P\\» h \\_i n < SUp dfh^hh + J £ cr a |g||| M/i _ p^f I . (5) 
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PROOF. Let fi h e A h . Then 

||/M||-i, 7 < llA*/i - + \\PfJ>h\\-^ n • (6) 

Using the definition of the norm on H~^(j), the fact that P is the orthogonal projection, Cauchy-Schwarz's 
inequality and the approximation properties of P (cf. [5]) it follows that 



\Hh- PfJ-hW-i , 7 = sup 



Ege 7 , (M/j - P^h,Oe 
SUP jr-jj 

sup 



< sup — * p| 

(E ge7 , lelll^ - ^WllU* {E?=i E, C73 I^IIC - PSWieY 

< sup ^ 

«6H^( 7 ) 



ll%7 



{Eee 7s C.|e|||^ - P^ll^} 5 {Ef=i ll?lll, 7 J ; 
< C sup — - 



< 



Cl J2Cs\e\\\^-P^\\h\ ■ (7) 



To bound the second term in ^ we start noting that using the continuous inf-sup condition (cf. |||) there 
exists (3 > such that 

P\\P» h \\_ hl < sup ~ ( ff'^ . (8) 

Next, since the pair Vh x satisfies a discrete inf-sup condition (cf. [B]) there exists a Fortin operator 
ith '■ Hq(Q) —> Vh, i.e., a continuous linear operator such that (/j, h ,v) 7 = (fi h , 7T/s(u)) 7 for all /i^ £ A,~. Then, 
using ([8]), the properties of 7r^, Cauchy-Schwarz's inequality, the approximation properties of P and the 
trace theorem it follows that 

3,1 S II / -(P^ h ,n h {v)) 7 
P \\PHh -i 7 < SUp — —- 

{Hh - P^h- / ^h{v)) 1 -(Vh,fth(v))-y 

< sup — — -—— + sup 



veH£(n) C\~*h(v)\i,n ven^n) C\^h(v)\i,n 

Eg 67 , Oft ~ PfJ-h, TTh{v) - P^h{v))e -U h , TT h (v)), 

< SUp —. 1-— j h SUp 



{E^aielll^-^II^I^E^II^WIIl^}^ 
< G sup — j rj f h C sup 



-{Hh-,Vh)r 



and the result follows. □ 



We next state the main stability result for For this, we introduce the following mesh-dependent 
norm on Wh- 

: = kli,n + P 2 \W\\l hl + E c M\»h - PfihWle ■ ( 9 ) 

Theorem 2. XTie bilinear form B satisfies 

sup rr-p > ^\\(uh,Xh)\\w h , (10) 

for all (uh,Xh) G Wft,. Hence, problem ^ is well-posed. 

Proof. Let (uh, A/j) € WV From the definition of £? it easily follows that 

B[(u h ,X h ),(u h ,X h )] = \u h \l n + E C a |e|||A fc -PA fc ||§ )g . (11) 

Next, from Lemma |lj there exists to/, G V/, such that ju^li^ = /3||A/j||_i 7 and 

/8 2 ||A/i||^i i7 -/3||A fc ||_i i7 j E C s \i\\\X h - PX h \\lA <-(\ h ,w h )t. 

Then, applying Cauchy-Schwarz's and Young's inequalities we obtain 

B[(uh, X h ), (uh + \ w h, A;,)] = \uh\\,n + E C s |e|||A/, - PX h \\l s + j[(Vii ft ,Vw A )sj - (Xh,w h )^\ 

> \Wh\ln + E C M\*h - PX h }\l s - ^/3 2 ||A,|| 2 _ i7 + ^l|Aft|| 2 _ i7 - \ E C.|e|||A„ - PX h f^ 
and the proof is finished noting that + §iOft, A/ l )||>v Jl < |||(u/j, A/ l )|jyy fe . □ 



4. Error analysis 

Supposing that A € £ 2 (7) we split the error into interpolation and discrete errors as follows (e u ,e\) := 
(u - u h , A - X h ) = (u - l h u, X - J h X) + (2 h u - u h , J h X - X h ) =: (77", ry A ) - (e£, e\), where l h stands for 
the Lagrange interpolation operator, and JhX € A/j is defined by JhX\ e := |e|~ 1 (A, l) e . As most of LPS-like 
methods, ^ introduces a consistency error. Using pj), ([3| and the definition of B, the following result is 
readily established. 

Lemma 3. Let us suppose that X E L 2 ^). Then, for all (vh, l^h) € Wft 

B[(e u ,e x ),(v h ,fi h )} = E C s \i\{X-PX^ h -P^ h )~ e . 

e£~t h 

Note that, though solution of ([I]) can be supposed in H 2 {w), we can not expect the same regularity for 
the solution u of ([2]), which only belongs to H s (£l), with | — e < s < 2 for any e > (see [6]). For the 
Lagrange multiplier A, it belongs to L 2 (j) in the worst case and to Tlf^H? (7^) in the best case. The main 
result of this section, namely the convergence of method ([3]), is stated next. 
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Theorem 4. Let us suppose that u € H s (Vt) - e < s < 2) and that A e n^ 1 iJ 5 (7 J ) (0 < S < \). Then 
there exists a constant C > 0, independent of h, such that: 



N 



(e u ,ex)\\w h < C \ h s - 1 \u\ s ^ + hi+ s [ £ ||A||| 



1i 



(12) 



Proof. The first step is to bound the discrete error. For this, let {wk^k) & W/, such that IKn^i^Hwh = 1 
and the maximum on Theorem [2] is attained. Then, using Lemma [3] and Cauchy-Schwarz's inequality we 
arrive at 

l\\(e h u ,e'm Wh < B[(e h u ,e h x ),(w h ,t h )} 
= - B[(e u ,e x ), (w h , t h )\ + B[(n u , r] X ), (w h) t h )\ 
= ~ E C s \e\{\ - P\t h - Pt h )- e + {Vr] u ,Vw h ) n - {^^h^ 

+ {th,V u )-y+J2 C s \e\(r) x -Pr) x ,t h -Pthh 
< brii.nKko + \\v X \\-iJw h \\ hl + ||*h||_i, 7 ||»7"||i l7 

+ V2 ( Cs\i\ [||A - PA||g, 8 + \\ V X Pv X \\h] ) ( E C M\*H PthWh 



Next, the fact that ||(w^,ift)|jvy h = 1, the trace Theorem and Poincare's inequality lead to 



\\WA)\\w*<C + ||r? A |U, 7 + 



E c >w 



|A-PA||g, 8 + h X -Pv X \\h 



(13) 



Using a standard interpolation result (cf. [5]), we have |ry M |i,n < Ch s 1 |tt| S) n. In addition, the approximation 
properties of P, its continuity and the approximation properties of lead to 

J2 IIA - PX\\le < C^||A||| 7j and £ ||^ - P^ ||* g < C7> 2<5 ||A||| 7 . . 
Finally, to bound the term ||?7 A ||_i 7 , we follow steps analogous to Q and use |e| < Ch to arrive at: 



(14) 



II \ a- Ml Sill (A - Jh\i~ JhOlj , ^il+i I ii \ II 

||A- JkA||_l. 7 = SUP r-r < Ch* +a [ 2J|A|| 



A' 



2 

5 .7i 



vi=i 



Hence, gathering (JJL3J) , ([15]) and ( |14[ ) we obtain 

h s - 1 \u\ St a + hi +s 



\\(e h u ,e h x )\\ Wh <C 



N 



En 



<5,7j 



The interpolation error \\(rj u , V )\\w h is bounded in the same way and we get 



|?fkn + /% A |U, 7 + 



E^l g IN A -^ A Ho,e| < 

V« 6 7s 



N 



h^lu^n + h^ 6 ( Ell A H 



2 

■5,73 



The error estimate then follows from the triangular inequality. 
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(15) 



□ 



5. A numerical experiment 




Figure 1: Convergence of the method |3|-(|4}: errors \u — «h|l,n (left) and ||A — || o,-y (right). 

In order to illustrate the above theoretical results, a numerical test has been performed. Problem 
([I]) has been solved using method (|3])-((4|- We chose uj = [0;1] 2 , VL = [— a; 1 + a] 2 (with a > 0), and 
f(x, y) = 2((x + a)(l + a — x) + (y + a)(l + a — y)) so that problem |T]) has an analytical solution u a (x, y) = 
(x + a)(l + a — x)(jj + a)(l + a — y). We set g = u a | 7 and use a = 0.5. A structured mesh Th of is built, from 
which the boundary meshes jh and 7^ were obtained automatically, with 7^ satisfying 3ft. < |e| < 6/1. The 
computations have been performed with Matlab™/ Octave. The errors appear to be fearly independent 
of the value of C s in the range 0.1 < C s < 1000, and then we have fixed C s = 0.1 in our experiments. 
Convergence results are displayed Figure]!] Note that for the errors \u — Uh\i,n and || A — A/jIIo^ the optimal 
convergence order O(h) is recovered, with a faster convergence rate for the latter (in our case A = 0, which 



helps to explain the faster convergence). This confirms our theoretical result (12). Without stabilization 
(C s = 0), a singular matrix is obtained if the condition 3h < \e\ is violated. This confirms the necessity 
of the geometric condition 3/i < |e| < Ch of without stabilization, and the interest of our stabilized 
formulation. 

Acknowledgements 

The second author acknowledges PEPII project SIM2E and LMB (UMR CNRS 6624) for partial funding. 



References 

[1] F. Brezzi, M. Fortin, A minimal stabilisation procedure for mixed finite element methods, Numer. Math. 89 (2001) 457-492. 
[2] E. Burman, Pressure projection stabilizations for Galerkin approximations of Stokes' and Darcy's problem, Numer. Methods 

Partial Differential Equations 24 (2008) 127-143. 
[3] E. Burman, P. Hansbo, Fictitious domain finite element methods using cut elements: I. A stabilized Lagrange multiplier 

method, Comput. Methods Appl. Mech. Engrg. 199 (2010) 2680-2686. 
[4] E. Burman, P. Hansbo, Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche's method, 

Appl. Numer. Math. (2011). To appear. 
[5] A. Ern, J.L. Guermond, Theory and practice of finite elements, volume 159 of Applied Mathematical Sciences, Springer- 

Verlag, New York, 2004. 

[6] V. Girault, R. Glowinski, Error analysis of a fictitious domain method applied to a Dirichlet problem, Japan J. Indust. 

Appl. Math. 12 (1995) 487-514. 
[7] R. Glowinski, T. Pan, J. Periaux, A fictitious domain method for Dirichlet problem and applications, Comput. Methods 

Appl. Mech. Engrg. Ill (1994) 283-303. 
[8] J. Haslinger, Y. Renard, A new fictitious domain approach inspired by the extended finite element method, SIAM J. Numer. 

Anal. 47 (2009) 1474-1499. 

[9] N. Moes, E. Bechet, M. Tourbier, Imposing Dirichlet boundary conditions in the extended finite element method, Int. J. 
Numer. Meth. Engrg. 67 (2006) 1641-1669. 



6 



